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Abstract 

We report a measurement of the charmless semileptonic decays -k^ /p~t^v and — > 

I (Pt^v, based on 253 fb~^ of data collected with the Belle detector at the KEKB e^e~ asym- 
metric collider. In this analysis, the accompanying B meson is reconstructed in the semileptonic 
mode B 0^*^11/, enabling detection of the signal modes with high purity. We measure the 
branching fractions tt'I+u) = (1.48 ± 0.20 ± 0.16 ± 0.04) x 10"'', B{B° p-i+u) = 

(2.07 ± 0.47 ± 0.25 ± 0.14) x lO"'', B{B+ ir^l+u) = (0.76 ± 0.13 ± 0.08 ± 0.04) x lO"^ and 
H(B+ p°£+u) = (1.39 dr 0.23 ± 0.17 ± 0.03) x 10"^, where the errors are statistical, experi- 
mental systematic, and systematic due to form-factor uncertainties, respectively. For each mode 
we also present the partial branching fractions in three q'^ intervals; q'^ < 8, 8 < q'^ < 16, and 
q^ > 16GcV^/c^. Based on these results, the magnitude of the Cabibbo-Kobayashi-Maskawa 
matrix element V^b is extracted. All of the presented results are preliminary. 

PACS numbers: 12.15.Hh, 12.38.Gc, 13.25.Hw 
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INTRODUCTION 



Exclusive B Xuii^ decays proceed dominantly via a 6 — > uW' tree process and can 
be used to determine \Vub\, one of the smallest and least known elements of the Cabibbo- 
Kobayashi-Maskawa matrix A major theoretical uncertainty stems from the specification 
of form-factors (FF) involved in the decays. The recent release of preliminary results from 
unquenched Lattice QCD(LQCD) calculations of FFs P, Isl] has made possible the model- 
independent determination of \Vub\- Since LQCD results are available only in the high 
region (> 16GeV^/c^), a clean measurement of the partial B — > Tciu branching fraction in 
the same high region is essential. 

There have been several measurements in the past by CLEO, BaBar and Belle for the 
B — > niu, piu, rjiu and uiu modes. The analyses in |3,H,M0,IEI3 utilize the method, origi- 
nally developed by CLEO, where the B decays are reconstructed by inferring the undetected 
neutrino momentum from missing energy and momentum ( 'V-reconstruction method" ) 0| . 
In the 5-factory era, we may straightforwardly improve the statistical precision by sim- 
ply applying the //-reconstruction method. However, the systematic uncertainty limits the 
measurement because of a poor signal-to-noise ratio. 

The Belle collaboration has previously presented measurements of B^ -k^ / p^l^v via 
D^*^^ decay tagging based on 140 fb~^ data The work demonstrated that the method 
allows one to detect the signal with high purity and is a promising approach to yield the best 
overall precision from the large data sample being accumulated by i?-factories over the next 
few years. The BaBar collaboration also presented measurements of exclusive B — > X„£z/ 
decays, both with D''*Hv decay tagging jll|, [l^ and with hadronic decay tagging [13^ . 

In this paper we present measurements of B^ t[~ / p~l^v and B^ —>■ vr^/p'^^+z/ decays 
using D^*^ih' decay tagging and based on 253 fb~^ data. In comparison with the previous 
measurement 3| , the analysis has been extended to higher integrated luminosity and to the 
additional modes B^ tt^/ p^i^u. The analysis technique is nearly identical to the previous 
one; we reconstruct the entire decay chain from the T(4S'), T(4S') BsigBtag, Bsig Tc/piu 
and Btag — ^ D^*^^ with several D^*^ sub-modes. The back-to-back correlation of the two 
B mesons in the T(4S') rest frame allows us to constrain the kinematics of the double 
semileptonic decay. Inclusion of charge conjugate decays is imphed throughout this paper. 



DATA SET AND EXPERIMENT 

The analysis is based on data recorded with the Belle detector at the asymmetric e"'"e~ 
collider KEKB operating at the center-of-mass (cm.) energy of the T(4S') resonance jl^ . 
The T(4S') data set used for this study corresponds to an integrated luminosity of 253 fb~^ 
and contains 275.2 x 10® BB events. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K^^ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 1151]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 152 million BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector, and a 
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small-cell inner drift chamber were used to record the remaining 123 million BB pairs |l6l|. 

A detailed Monte Carlo (MC) simulation, which fully describes the detector geometry and 
response and is based on GEANT jl7|J, is applied to estimate the signal detection efficiency 
and to study the background. To examine the FF dependence, MC samples for the B — > Tiiu 
signal decays are generated with different form-factor models; a quark model (ISGW II jlS^). 
light cone sum rules (LCSR; Ball'Ol for vr^z/ [l^ and Ball'QS for piu ^2d\) and quenched 
lattice QCD (UKQCD [2l|). A relativistic quark model (Melikhov |22|]) is also used for 
piu. To model the cross-feed from other B —>■ Xud-v decays, MC samples are generated with 
the ISGW II model for the resonant components and the DeFazio-Neubert model [l^l for 
non- resonant component (vrfz/ and plv components are excluded in this sample). To model 
the BB and continuum backgrounds, large generic BB (based on Evtgen j23]) and qq Monte 
Carlo samples are used. 



EVENT RECONSTRUCTION AND SELECTION 

Charged particle tracks are reconstructed from hits in the SVD and CDC. They are re- 
quired to satisfy track quality cuts based on their impact parameters relative to the measured 
profile of the interaction point (IP profile) of the two beams. Charged kaons are identified 
by combining information on ionization loss (dE/dx) in the CDC, Cherenkov light yields 
in the ACC and time-of-fiight measured by the TOF system. For the nominal requirement, 
the kaon identification efficiency is approximately 88% and the rate for misidentification of 
pions as kaons is about 8%. Hadron tracks that are not identified as kaons are treated as 
pions. Tracks satisfying the lepton identification criteria, as described later, are removed 
from consideration. 

Neutral pions are reconstructed using 7 pairs with an invariant mass between 117 and 
150 MeV/c^. Each 7 is required to have a minimum energy deposit of 50 MeV. Kg mesons are 
reconstructed using pairs of charged tracks that have an invariant mass within ±12MeV/c^ 
of the known Kg mass. 

Electron identification is based on a combination of dE/dx in CDC, the response of ACC, 
shower shape in ECL and the ratio of energy deposit in ECL to the momentum measured by 
the tracking system. Muon identification by KLM is performed by resistive plate counters 
interleaved in the iron yoke. The lepton identification efficiencies are estimated to be about 
90% for both electrons and muons in the momentum region above 1.2 GeV/c. where leptons 
from prompt B decays dominate. The hadron misidentification rate is measured using 
reconstructed Kg tt+tt^ and found to be less than 0.2% for electrons and 1.5% for muons 
in the same momentum region. 

For the reconstruction of Btag ~^ D^*^£i', the lepton candidate is required to have the 
correct sign charge with respect to the D meson fiavor and a laboratory momentum greater 
than l.OGeV/c (p^"^ > 1.0 GeV/c). The D meson candidates are reconstructed by us- 
ing 7 decay modes of - — >• K'tt^ti^^ K'-k^-k^-k^, K'gTi^ , K^ti^-k^^ ii'^Tr+Tr+Tr", 
K+K% K+K-71+- and 10 decay modes of - ^ K-n+, K-n+7i^, K-tt+ti+ti-, /sT^tt^, 
K^gTi+Ti-, K^gii+TX-TX^, i^-7r+7r+7r-7r°, K+K", K^gK+K-, KgK-7r+. The candidates are 
required to have an invariant mass rriD within ±2.5cr (a is a standard deviation) of the 
nominal D mass, where the mass resolution a is dependent on the decay mode. D* mesons 
are reconstructed by combining the D meson candidate and a charged or neutral pion, 
D*~^ D^Tr"*", Z)"'"7r° and D*'^ D^tt^. Each D* candidate is required to have a mass 
difference Am = — nif) within ±2.5(t of the nominal values. 
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For the reconstruction of Bsig Xuijy, the lepton candidate is required to have the right 
sign charge with respect to the X„ system and p^"^ > 0.8 GeV/c. The Xu system may 
consist of one pion or two pions ( Nt,+ = 1 or iV^+ = iV^o = 1 for a tag and N^^o = 1 
or Nt,+ = Nt^~ = 1 for a B~ tag). The event is required to have no additional charged 
tracks or 7r° candidates. We also require that the residual energy from neutral clusters be 
less than 0.15 GeV {Eneut < 0.15 GeV). The two leptons on the tag and the signal sides are 
required to have opposite charge. The loss of signal due to B^ — B^ mixing is estimated by 
MC simulation. 

We then impose a constraint based on the kinematics of the double semileptonic decay 
in the T(4S') rest frame. In the semileptonic decay on each side, -Bi(2) Yi{2)^ (^i = D'^*H 
and Y2 = Xui), the angle between the -81(2) meson and the detected Yi(2) system Ob^^) 
is calculated from the relation, pi = {ps — Vy)"^ = and the known Pg (the absolute 
momentum of the mother B meson). This means that the -81(2) direction is constrained 
on the surface of a cone defined with the angle Ob^^) around the direction of the Y'i(2) 
system, as shown graphically in Figure ^ Then the back-to-back relation of the two B 
meson directions implies that the real B direction is on the cross lines of the two cones 
when one of the B system is spatially inverted. Denoting 612 the angle between the D*l 
and the X^f systems, the B directional vector ub = {xB,yB, zb) is given by, zb = cosObi, 
Ub = {cos9b2 — cos6'B2'^os6'i2)/sin6'i2, and 



xb = ±1/1 . (cos^Obi + cos^6b2 — 2cos6bj^cos9b2Cos9i2) (1) 

y sm6'i2 

with the coordinate definition in Figure If the hypothesis of the double semileptonic 
decay is correct and all the decay products are detected except for the two neutrinos, x% 
must range from to 1. Events passing a rather loose cut x"^ > —2.0 are used for signal 
extraction at a later stage of the analysis. 



nB = ( XB, yB, ZB ) 




X 



^PY2 



FIG. 1: Kinematics of the double semileptonic decay. 



Since Eq.^has two solutions and the direction of the B meson is not uniquely determined, 
we calculate, as = [E;^^^- E*^y - p\^^ , using the beam energy {El^^^), energy {E*^^) 
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and momentum (Px^) of the Xu system and neglecting the momentum of the B meson in 
the cm. system. The signal Monte Carlo simulation predicts the average resolution to 
be approximately 0.75GeV^/c^. 

According to Monte Carlo simulations, the major backgrounds originate with B —* X^iu 
and non-signal B X^^v decays, where some particles escape detection. There are sizable 
contributions from cross talk between the 5° and B^ tags. The contribution from qq 
processes is found to be negligible. 

For events selected as described above, the signal MC simulation indicates that the total 
detection efficiency (etotai), on the average of the electron and muon channels, is 2.00 x 10~^ 
for n-£+u and 7.75 x 10"^ for p~£+z/, 1.49 x lO'^ for 7r°£+z/ and 1.76 x lO'^ for p^i+u 
assuming the LCSR model. Here, etotai is defined with respect to the number of produced 
BB pairs where one B decays into the signal mode, and includes the loss of signal due to 
B^ — B^ mixing. Because of the relaxed lepton momentum cut (> 0.8 GeV/c), the variation 
of efficiency with different FF models is relatively small. 

The validity of the method to reconstruct the double semileptonic decay is checked by 
reconstructing the signal as B^^^ — > D*~i^i' followed by D*~ — > D^tt~,D^ K~^n~, with 
the same requirement on the tagging side. Figure |21-a) shows the obtained Mk-k-k distribution 
and its MC expectation. With the 253 fb~^ data sample, 224.7 ± 15.4 decays are clearly 
identified, while the expectation based on the MC efficiency and the product branching 
fraction B{B^ D*-i+u,D*- DH-,!)^ K+n-) = (1.38 ± 0.06) x 10"^ is 
224.5 ± 9.5 events. Their ratio R = 1.001 ± 0.093 is consistent with unity within the error, 
where 8.3% originates from the statistics of the detected D*~l^v decays and 4.9% from 
the error on B{B^ D*^l^v). Figure |21-b) shows the comparison of the reconstructed x\ 
distribution using the same sample to the MC simulation. The agreement between the data 
and MC demonstrates the validity of the present measurement. 




FIG. 2: Reconstructed M{KTnr) distribution(left) and distribution(right) for the D* i^i/ 
calibration decay. Points with error bars are data and the histogram is the signal MC. 
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EXTRACTION OF BRANCHING FRACTIONS 

The -K^ / p^l^v and -n^ j p^t^v signals are extracted by fitting the obtained 

two-dimensional distribution in (x^,Mx), where Mx is the invariant mass of the X„ sys- 
tem. Here, the fit components are the four signal modes, the other B^ X^i'^v and 
B^ X^t^v backgrounds and the BB background (seven components in total). The PDF 
(probability distribution function) for each fit component is determined from MC simula- 
tion. The -nj pt^v signal events exhibit characteristic behavior in both of their x\ and Mx 
distributions; other B X^t^v events exhibit a weak peaking structure in x\ but a broad 
distribution in Mx] the BB background has a relatively flat distribution in x\ and a broad 
structure in Mx- The PDFs in (x^, Mx) for each of the seven fit components are prepared 
for both B^ and B~ tag candidates. We then fit the two (a;^,Mx) distributions simulta- 
neously, constraining contributions from the cross talk between B^ and B~ tags as well as 
the cross-feed between vrf+z/ and pt^v. The fitting is constrained so that the sum of the 
deduced branching fractions for B iri^u, B pC^v and B other X„£+z/ is equal to 
the total inclusive branching fraction B{B XJv) = (0.25 ± 0.06)% f^. 

Figure El presents the projections on Mx and of the fitting result for data in the entire 

region. The extracted yields for the signal components are N{7T~i~^i') = 150.7 ± 19.7, 
N{p-i+u) = 88.6 ± 18.9, N{n^i+u) = 71.7 ± 11.5 and N{pH+u) = 143.3 ± 25.4, with the 
LCSR model used for the four signal PDF. 

Figure m shows projections of the data, separated into three bins, < 8 GeV^/c^ and 
8 < g^ < 16 GeV^/c^, g^ > 16 GeV^/c^. Here the normalization of the other B — > Xuiu 
and the BB background components are fixed to those obtained in the above fitting for 
the entire g^ region. Table Q] summarizes the extracted branching fractions with different 
FF-models for each signal mode, respectively. The results are unfolded using an efficiency 
matrix that relates the true and reconstructed g^ for the three g^ intervals. We calculate 
the total branching fraction by taking sum of the partial branching fractions in the three g^ 
intervals. 

SYSTEMATIC ERRORS 

Tables |H] and IIIII summarize the experimental systematic errors on the branching frac- 
tions. The experimental systematic errors can be categorized as originating from uncertain- 
ties in the signal reconstruction efficiency, the background estimation, and the normalization. 
The total experimental systematic error is the quadratic sum of all individual ones. We also 
consider the systematic error due to the dependence of the obtained branching fractions on 
the FF model. 

The effect from the uncertainty on the signal reconstruction efficiency is evaluated based 
on the efficiency calibration with the Bgig D*^i~^i> sample, discussed above. The error is 
taken to be that on the ratio of observed to expected number of the calibration signals (9.3%). 
This gives the largest contribution to the systematic error. Note this error is dominated 
by the statistics of the calibration signals, as explained above. Therefore, accumulation 
of additional integrated luminosity in the future will help to reduce this uncertainty. We 
further include residual errors for the reconstruction of the signal side; 1% and 2% for the 
detection of each charged and neutral pion, respectively, and 2% for the kaon veto and 2.1% 
for the lepton selection. 

The systematic error due to the uncertainty on the inclusive branching fraction B{B 
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FIG. 3: Projected Mx distribution(left) and distributions for the mass region of vr {Mx < 
0.18GeV/c^, middle) and p (0.6 < Mx < l.OGeV/c^, right) in all region; points are data. 
Histogram components are TT^i^u (red narrow 135° hatch), p~i~^u (red wide 45° hatch), other 
Xui~^i^ from (red cross-hatch) and TT^i^u (blue narrow 45° hatch), p^i'^u (blue wide 135° 
hatch), other Xui'^i' from (blue cross-hatch) and BB background (green border). 



Xuiiy), which is used to constrain X~i~^i' background, is estimated by varying this parameter 
within ±la of the error. The uncertainty on the BB background shape in our selection cut 
{Nt,+ = 1 or = N^o = 1 for a 5° tag and A^^o = 1 or A^^+ = N^- = 1 for a B~ tag) is 
studied in the simulation by randomly removing charged tracks and 7r° according to the error 
in detection efficiency (1% for a charged track, 2% for vr"), and also by reassigning identified 
charged kaons as pions according to the uncertainty in the kaon identification efficiency 
(2%). The resultant changes in the extracted branching fractions are assigned as systematic 
errors. We have seen significant change for pi'^u in the high region (g^ > 16 GeV^/c^). 
We also varv the fraction of 5 — > D**iv decays in the BB background MC by the error 
quoted in [25| to test the B Xc^l' model dependence in the BB background shape. To 
assess the uncertainty due to the production rate of we vary the production rate by the 
error of the inclusive branching fraction oi B ^ X quoted in 2^ . 

As for the normalization, we consider the uncertainty in the number of B^B^ and B^B~ 
pairs; the ratio of B^B~ to B^B^ pairs {f^/ fo, 2.5%), the mixing parameter {xd, 1.0%) and 
the measured number of BB pairs {N^g, 1.1%). 

The dependence of the extracted branching fractions on the FF model has been studied 
by repeating the above fitting procedure with various FF models for the signal mode and 
also for the cross- feed mode (vr^z/ ^ piu) . We consider the models listed in Table B For 
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FIG. 4: Projected Mx distribution in each region; Dot is data. Histogram components are 
■n^l^v (red narrow 135° hatch), p~i^v (red wide 45° hatch), other Xul^v from (red cross- 
hatch) and TT^l'^i' (blue narrow 45° hatch), p^l'^u (blue wide 135° hatch), other X^f.^i' from 
(blue cross-hatch) and BB background (green border). 



the extracted B{B -k^ the standard deviation among the models is < 1.7(4.4)% 
for T[l^v and < 1.9(1.4)% for pi'^u. For B{B —>■ p"(p°)£"'"z/), the standard deviation is 
< 6.2(1.8)% for pi'^u and < 2.5(0.1)% for ni'^i'. The assigned total error due to the FF 
model dependence is the quadratic sum of the maximum variations with the signal and 
cross-feed FF models. 



RESULTS 



Table IIVI summarizes our measurements of the total and partial branching fractions for 
the four signal modes. Each branching fraction is obtained by taking the simple average of 
the values obtained from the FF models shown in Table lU The errors shown in the table 
are statistical, experimental systematic, and systematic due to form-factor uncertainties. 
The obtained branching fractions for — »• 7r~_lp~ u are consistent with the existing 
measurements by CLEO 0| and BaBar 0, 0, EH, E^, within the measurement uncertainties. 
The overall uncertainty on our result for 5° — >■ 7r~£"'"z/ is 18%, comparable to that on the 
measurement published by CLEO ^ which is based on //-reconstruction. 

Figure presents the measured distributions for each signal mode, overlaid with the 
best fits of FF shapes to the data. To be self-consistent, the shape of a particular FF model is 
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TABLE I: Extracted branching fractions for each signal mode with different FF models in units 
of 10~^; the total branching fraction and the partial branching fractions in three intervals. 
X^/n and the probability of shows the quality of the fit of the FF shape to the extracted 
distribution. 



Mode 



Model 



BairOl 
ISGW II 
UKQCD 
Average 



BairOS 
ISGW II 
UKQCD 
Melikhov 
Average 



BairOl 
ISGW II 
UKQCD 
Average 



Bair98 
ISGW II 
UKQCD 
Melikhov 
Average 



total 



B 



<8 



8-16 



B 



>16 



Prob. 



1.48 ±0.20 
1.45 ±0.20 
1.50 ±0.20 
1.48 ±0.20 



0.56 ±0.11 
0.53 ±0.11 
0.55 ±0.11 
0.55 ±0.11 



0.51 ±0.11 
0.53 ±0.11 
0.55 ±0.12 
0.53 ±0.12 



0.41 ±0.12 
0.39 ±0.11 
0.41 ±0.11 
0.40 ±0.12 



2.7/3-1 
4.2/3-1 
0.0/3-1 



0.26 
0.13 
0.98 



2.08 ±0.47 

2.06 ±0.47 
1.97 ±0.44 
2.20 ±0.49 

2.07 ±0.47 



0.51 ±0.22 
0.50 ±0.20 
0.48 ±0.21 
0.50 ±0.21 
0.50 ±0.21 



1.07 ±0.29 

1.00 ±0.28 

1.01 ±0.28 

1.02 ±0.29 
1.02 ±0.29 



0.49 ±0.27 
0.59 ±0.30 
0.48 ± 0.27 
0.68 ± 0.34 
0.56 ± 0.30 



2.1/3-1 
0.6/3-1 
1.0/3-1 
3.3/3-1 



0.35 
0.73 
0.59 
0.19 



0.74 ±0.13 
0.76 ±0.13 
0.80 ±0.14 
0.76 ±0.13 



0.21 ±0.06 
0.20 ±0.06 
0.22 ±0.07 
0.21 ±0.06 



0.32 ±0.08 
0.33 ± 0.09 
0.35 ±0.09 
0.33 ± 0.09 



0.21 ±0.08 
0.22 ±0.08 
0.23 ±0.08 
0.22 ±0.08 



2.0/3-1 
6.6/3-1 
1.1/3-1 



0.36 
0.04 
0.58 



1.40 ±0.23 

1.38 ±0.23 
1.36 ±0.23 
1.40 ±0.23 

1.39 ±0.23 



0.50 ±0.12 
0.48 ±0.12 
0.48 ±0.12 
0.51 ±0.12 
0.49 ±0.12 



0.57 ±0.15 
0.56 ±0.15 
0.58 ±0.15 
0.57 ±0.15 
0.57 ±0.15 



0.33 ±0.12 
0.34 ±0.12 
0.31 ±0.12 
0.32 ±0.12 
0.32 ±0.12 



1.5/3-1 
3.2/3-1 
1.4/3-1 
1.9/3-1 



0.48 
0.21 
0.49 
0.39 



fit to the distribution extracted with the same FF model. The quality of the fit in terms 
of the probability of x^? shown in Table U may provide one way to discriminate 

among the models. For both —>■ vr^f+z/ and —>■ n^i^i' the agreement with the ISGW 
II model is marginal, as indicated by other measurements p, 0|- 

In this work, the 5° ^ 7r-£+i//5+ n^i+u and 5° ^ p-i+u/B+ p^i+u signals 
are extracted separately, which allows us to test the isospin relations. /^From the obtained 
branching fractions and the B meson lifetime in , the ratios of decay rates are found to 
be, 



F(50 ^ n-i+u) 
r{B+ 7rO£+z/) 

T{B+ pH+u) 



(2.08 ±0.45 ±0.17), (2) 
(1.59 ±0.44 ±0.17), (3) 



where the first and second errors are statistical and systematic errors, respectively. The 
both ratios are found to be consistent with the isospin relations; F(i?° 7r^(p^)£+z/) = 
2F(S+ ^ 7rO(pO)£+zy). 

The obtained branching fractions in Table IIVI can be used to extract \Vub\ using the 
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TABLE II: Summary of systematic errors(%) for tt / p l^v. 



Source 


g2 < 8 


8 - 16 > 16 


< 16 


all 


g2 < 8 


p-l^u 
8 - 16 > 16 


< 16 


all 


Tracking efficiency 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


7r° reconstruction 


- 






- 


- 


2 


2 


2 


2 


2 


Lepton identification 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


Kaon identification 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


D*iv calibration 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


Br{Xuf.v) in the fitting 


0.2 


2.0 


2.1 


1.1 


1.4 


0.8 


4.2 


13.1 


1.1 


5.4 


BB background shape 


1.5 


1.7 


3.5 


1.3 


1.7 


2.5 


1.8 


9.0 


1.3 


2.6 


Br{D**lu) 


1.2 


1.4 


0.1 


1.4 


0.9 


0.5 


0.1 


1.3 


1.4 


0.4 


Kj^ production rate 


0.2 


0.3 


0.4 


0.3 


0.3 


1.0 


0.8 


2.2 


0.3 


1.2 




1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


f+/fo 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


Xd 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


exp. total 


10.4 


10.7 


11.0 


10.5 


10.5 


10.8 


11.4 


19.2 


10.5 


12.1 


F.F for signal 

F.F for cross-feed 


2.6 

2.2 


4.0 

3.9 


3.1 

5.1 


1.9 
2.8 


1.7 
1.9 


3.3 

2.8 


5.4 

2.0 


21.4 

6.6 


1.9 

2.8 


6.2 

2.5 


F.F total 


3.4 


5.6 


6.0 


3.3 


2.5 


4.3 


5.8 


22.4 


4.8 


6.7 



TABLE III: Summary of systematic errors(%) for B~^ 'k^ j fPl'^v. 







7T 










p^i+f 






Source 


g2 < 8 


8 - 16 > 16 


< 16 


all 


g2 < 8 8 


- 16 


> 16 


< 16 


all 


Tracking efficiency 












2 


2 


2 


2 


2 


7r° reconstruction 


2 


2 


2 


2 


2 












Lepton identification 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


Kaon identification 












4 


4 


4 


4 


4 


D*iv calibration 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


Br{X.,J-i^) in the fitting 


0.4 


3.7 


0.5 


1.3 


1.6 


2.3 


2.5 


10.8 


2.5 


4.4 


BB background shape 


2.2 


5.6 


2.7 


2.3 


3.8 


5.9 


4.3 


12.9 


1.7 


1.9 


Br{D**iiy) 


1.3 


0.7 


0.8 


0.5 


0.9 


0.2 


1.7 


3.1 


0.9 


1.4 


production rate 


0.4 


1.1 


0.7 


0.5 


0.8 


0.5 


1.1 


1.8 


0.2 


0.3 


^BB 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


f+/fo 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


exp. total 


10.5 


12.2 


10.5 


10.5 


11.0 


12.6 


12.1 


20.5 


11.3 


12.0 


F.F for signal 


5.1 


4.5 


4.4 


4.7 


4.4 


2.8 


1.3 


6.1 


1.7 


1.8 


F.F for cross-feed 


0.1 


2.1 


2.6 


1.6 


1.4 


1.0 


1.3 


4.4 


1.1 


0.1 


F.F total 


5.1 


5.0 


5.1 


5.0 


4.6 


2.9 


1.8 


7.5 


2.0 


1.8 
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TABLE IV: Summary of the obtained branching fractions. The errors are statistical, experimental 
systematic, and systematic due to form-factor uncertainties. 



Modes ref 


;ion (GeVVc^) 


Branching fraction (xlO ^) 






Total 


1.48 ±0.20 ±0.16 ±0.04 






> Id 


0.4U ± U.iz ± 0.04 ± 0.02 






< 16 


1.08 ±0.16 ±0.11 ±0.04 


B+ 




Total 


0.76 ±0.13 ±0.08 ±0.04 






> 16 


0.22 ±0.08 ±0.02 ±0.01 






< 16 


0.54 ±0.11 ±0.06 ±0.03 


SO 




Total 


2.07 ±0.47 ±0.25 ±0.14 


B+ 




Total 


1.39 ±0.23 ±0.17 ±0.02 




p'Tv 




HISGW2 
BallOl 
■ UKQCD 
B Melikhov 



a 



q^(GeVT 




q'(GeVT 



q^(GeVT 



FIG. 5: Extracted (f' distribution. Data points are shown for different FF models used to estimate 
the detection efficiency. Lines are for the best fit of the FF shapes to the obtained distribution. 
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relation, 

ivj . ME^, (4) 

\ -L thy Tb 

where Tthy is the form-factor normahzation, predicted from theories. We hst the partial 
branching fractions for B ttI^v decays in the region above 16 GeV^/c^, where the 
LQCD calculations are the most sensitive. The table provides also the results in the region 
below 16 GeVVc^, so that one can deduce \Vub\ based on other approaches such as the LCSR 
calculations isl l. 



In this paper, we quote \Vub\ based on the t^I^v data in the high region and the form 
factor predicted by recent unquenched LQCD calculations. Their predicted Vthy for the 

> 16 GeVVc^ region are fthy = 1-31 ± 0.33 ps~^ (HPQCD) |^ and 1.83 ± 0.50 ps"^ 
(FNAL) 0. We use tbo = 1.532 ± 0.009 ps and tb+ = 1.638 ± 0.011 ps and we 

assume the isospin relations for Tthy between tt^ and B^ ir^ transitions. Table IVl 

summarize the results, where the first and second errors are the experimental statistical and 
systematic errors, respectively. The third error is based on the error on Tthy quoted by the 
LQCD authors. These theoretical errors are asymmetric because we assign them by taking 
the variation in \ Vub\ when Tthy is varied by the quoted errors. By taking the average of the 
results obtained from B^ n^i^u and B^ ir^i^u data, we obtain 

iVubrig^iT''^"^ = (4-50 ± 0.52 ± 0.27^°:!°) x 10-3(HPQCD), (5) 
iKfeltSr"'^' = (3-81 ± 0.44 ± 0.23^^3) x 10~=^(FNAL). (6) 



TABLE V: Summary of \Vub\ obtained from the B ni'^i' data in the > 16GeV^/c'^ region. 
The first and second errors are experimental statistical and systematic errors, respectively. The 
third error stems from the error on T^hy quoted by the LQCD authors. 



Theory Tthyjvs-^) Mode |Kb|(xlO-^) 

HPQCD 1.31 ±0.33 n-i+u 4.49 ± 0.65 ± 0.28t^:J^ 

TT^£+u 4.53 ± 0.81 ± O.lQtoll 
■K^l+v + TT^l+u 4.50 ± 0.52 ± 0.27j:[J;I^ 
FNAL 1.83 ±0.50 -k-I+u 3.80 ± 0.55 ± 0.24t°j| 

vrO^+i/ 3.83 ± 0.69 ± 0.22l° 
■K-l+v + -K^l+u 3.81 ± 0.44 ± 0.23l|]:^^ 



These values are in agreement with those from the inclusive B Xuijy decays [2^ . The 
experimental precision in the above \ Vub\ determination is 13%, and it is currently dominated 
by the statistical error of 12%. By accumulating more luminosity, a measurement with errors 
below 10% is feasible. With improvement in unquenched LQCD calculations, the present 
measurement will provide a more accurate determination of \Vub\- 
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